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Abstract 

We re-analyze a possible ambiguity in the application of dimensional regulariza- 
tion to Einstein-Gauss-Bonnet gravity, arising from the way one treats the Gauss- 
Bonnet term T. It is demonstrated that the addition of such a term to the action 
gives rise to a non-minimal graviton wave operator, but does not produce new on 
shell divergences at one loop order in d = 4. However, from a d-dimensional per- 
spective the Gauss-Bonnet term is shown to generate new divergences in the form 
of the six-dimensional Euler density. The conjecture that one would next produce 
the eight-dimensional Euler term is shown to be false. 

1 Introduction 

The natural generalization of the Einstein-Hilbert action to a spacetime of dimension 
greater than four is provided by the Lovelock action [21 E]; see e.g. [1] for a review. It 
consists of the sum of all dimensionally continued Euler densities, each term coming with 
its own coupling constant. The Lovelock action gives rise to a symmetric conserved field 
equation which contains derivatives of the metric no higher than second order and is quasi- 
linear in these second derivatives. This leads to a well-defined classical Cauchy problem [51 
E], and to unitarity at the quantum level, unlike generic higher derivative actions [71IH]. 
Unitarity is also the reason why one expects the Lovelock action to make its appearance 
in the low-energy limit of string theory [9l llOj . More recently, numerous papers have 
appeared in which this same action appears in connection with brane scenarios in string 
theory, see e.g. [^ 1111 . This partly motivated our examination of the ultraviolet behavior 
of the Lovelock action. 

It is well known that a perturbative treatment of Einstein gravity gives rise to a one 
loop finite S-matrix in four dimensions. The background field method and dimensional 
regularization predict one loop divergences of curvature-squared form [12'. Although 
these do not have the functional form of the classical action, the vacuum field equa- 
tions combined with the Gauss-Bonnet identity would imply that such divergences do 
not survive. However, this reasoning was criticized in a little known paper by Capper 
and Kimber They pointed out that dimensional regularization requires one to work 
consistently in d dimensions and, although the field equation can trivially be continued. 
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there does not exist a d-dimensional Gauss-Bonnet identity. In short, the d — > 4 hmit of 
the expression 



is ill defined. This motivated Capper and Kimber to add a Riemann-squared term to the 
classical action. They then demonstrated that this extra term has no effect on tree-level 
graviton-graviton scattering. This was shown to be not due to explicit factors of d — 4, 
as none appear, but rather caused by the imposition of on shell conditions and setting 
d equal to four for all external legs. Although no one loop calculation was performed, 
an appeal to consistency with other regularization schemes led them to conclude that 
the one loop S-matrix of standard gravity is finite after all. The re-examination in the 
present paper answers the question of Capper and Kimber for the divergent parts of 
the effective action at one loop level. No changes are found, so we confirm what one 
would naively have expected from adding such a, d — A topological term. However, in the 
process we discover that from the c?-dimensional point of view advocated by Capper and 
Kimber, adding the Gauss-Bonnet term to the classical action next requires the d = 6 
Euler density as a counter term. Hence, it would seem that the rcnormalization process 
induces the full Lovelock action. 

Besides the well known general difhculty of any calculation in quantum gravity, a 
further reason why the problem was not yet studied beyond tree level is that one is then 
faced with a non- minimal wave operator. If the leading part of a wave operator equals 
the Laplacian, one speaks of a minimal operator. In gauge theories this form can usually 
be arranged by a judicious gauge choice. Minimal wave operators allow the use of the 
powerful Schwinger-DeWitt method and convenient background field algorithms are then 
available. We will show that the addition of the Gauss-Bonnet term inevitably leads to 
a non-minimal wave operator for the graviton. An elegant extension of the Schwinger- 
DeWitt method to the non-minimal case was given in [13j but we did not rely on this 
work. Since the offending term in the wave operator turns out to be of first order in the 
curvature, we may treat it perturbatively and thus return to the minimal setting. Using 
the fully covariant Schwinger-DeWitt method, we reproduce and extend to curved space 
an earlier flat space algorithm [14] . There, 't Hooft's background field algorithm [15] was 
generalized to second order in the non-minimal part of the wave operator (see also the 
recent systematic work of [TH ]T7] ) . 

In related work, Berredo-Peixoto and Shapiro recently demonstrated that no new one 
loop divergences are generated upon adding the Gauss-Bonnet term to either confor- 
mal [IH] or general [IH] higher derivative gravity. Such theories are known to be renor- 
malizable, though not unitary. The authors of [18|, 119] employed dimensional regular- 
ization and the Schwinger-DeWitt method for their involved calculations. Interestingly, 
the Gauss-Bonnet term was shown to actually affect the d-dimensional renormalization 
group equations and new non-trivial fixed points were found. However, these papers do 
not answer the original issue of [T] : Although the classical action in [T^] includes Einstein- 
Hilbert and cosmological terms, a higher derivative gauge fixing term was chosen which 
does not allow one to regain the special case of two-derivative gravity. Hence, the analysis 
in [12] excludes the case of Lovelock gravity we are interested in. 
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2 Perturbative expansion of the Lovelock action 



In a d dimensional Riemann space, the Lovelock action is given by 
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For compactness, the Riemann curvature tensor R^^'^" has been written here as i?^^ 
and dv = d'^x^. In d = 2k, the integrand of Sk is proportional to the Euler density. 
Note the total antisymmetry on 2k indices. For given dimension d, the Schouten identity 
implies that the maximum value of k is the integer part of d/2, denoted here as [d/2]. 
The Afc are arbitrary coupling constants. In four dimensions this action equals 
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dv[Xo + R + \2{R'"''"' R^..pa - ^R'^'Rt.. + R^) ] 



(2) 



One recognizes here the cosmological constant, Einstein-Hilbert term, and four-dimensional 
Euler density, also known as Gauss-Bonnet term. The latter is only locally a total deriva- 
tive and cannot be written as the divergence of a vector field. In higher dimensions, 
further dimensionally continued Euler densities appear. For later reference we note that 
in a Ricci flat space 5*3 reduces to (C represents the Weyl tensor) 
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dv {C\ 
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(3) 



Here, the labels W and M stand for Wheel and Mobius, respectively. Using the graphical 
notation of the Weyl-tensor shown in figure [1] one gets a simple graphical representation 
of C'^ and (see figure [U 



y jj, 



> — < 




Figure 1: Graphical notation for the Weyl-tensor and its symmetries. 
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Figure 2: Graphical representation of the Wheel and Mobius graphs 



In order to set the stage for our one loop calculation, we now make the usual background- 
quantum splitting via the replacement g^i, — s- g^^ + nh^^, where = 32TrG. The action 
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5 of Eq (IT|) then splits into k" ^5^"), the superscript n indicating the order in the 
quantum field h. We find 

. [(d-l)/2] 
k=Q 

Note that this expression is totally antisymmetric on 2fc + 1 indices, so on one more index 
than in Eq ((T|). The field equation can be read off directly from Eq E.g. in d = 6 

Upon writing out the indicated antisymmetrizations, one recognizes here the Einstein 
and Bach-Lanczos tensors in the second and third term, respectively. Observe that the 
field equation is covariantly conserved: Applying and using the Bianchi identity, each 
term vanishes separately, due to the invariance under general coordinate transformations. 
The contracted field equation reads 

dXo + id- 2)R + (d - A)\2{R^''P^ R^.pa - ^R'^'Rt.. + = (6) 
At second order in the quantum fields, we find 

„ [{d-l)/2\ 
•' k=0 

+ 'Vl -^A'2M3 ^^2(0-2^2(0-1 VV2fc;A'2fc + l 2-^^21= a '^tJ.2k + l)\ ^ ' 

where we used semi-colon to denote covariant differentiation. Explicitly in c? = 6 

S^-^^ = J dv^[\o{h^:h^;^-^h^) 

+3((/i>[r - khhl^)R''4^ + '^h^UhVf) + h^HR^Jh-^p) (8) 
+30A2((/i;^/i[r - ^hhl^)R-CPRaP + ^h^HR'C'phV-P + 2h^llRlPR7c}K)] 

We emphasize that the form of S''-^) is completely determined by the requirement of 
total antisymmetry in its indices; our calculations are only needed to find the numerical 
coefficients in Eq ([5]). We also note that the covariant derivatives in the penultimate term 
of Eq ([HI) can be moved via partial integration from one field h to another by grace of the 
Bianchi identity and the indicated total antisymmetry. The last line of Eq ^ shows that 
in a general curved background the Gauss-Bonnet term docs contribute to the quadratic 
terms and hence to the wave operator. Only in special backgrounds, e.g. in a flat space, 
does its contribution vanish. Eq ([5]) can be written as J h/S.h where the wave operator 
takes the schematic form 

A - ^ Afci?'=-iVV + XkR" (9) 

k=l k=Q 

Here, R represents the Riemann tensor or any of its contractions. In the next section 
we will show that the leading term of the first sum can be gauge fixed to the Laplacian 
but we will also demonstrate that the Gauss-Bonnet term unavoidably makes the full 
wave operator non-minimal. In maximally symmetric spaces the wave operator reduces 
to the minimal form (1 + cX2R)0 with R the Ricci scalar and c a number. This fact was 
exploited in [50] ■ We will not assume such a special background. 
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3 Gauge fixing the Einstein- Gauss-Bonnet action 



Before we can do quantum calculations, we need to fix the gauge. We choose the usual 
background covariant harmonic, or DeWitt-Feynman, gauge. We do so by adding the 
following gauge fixing term to the classical action in Eq ^ 

Sfr. = Jdvg^'^F^F, , - Z^;., - i/i.^ (10) 

Note that this will gauge fix the Einstein-Hilbert term but not the other terms. This 
suffices to define the propagator but does not affect the contributions coming from the 
Gauss-Bonnet term. Hence, we will wind up with a non-minimal wave operator. One 
might think that a more clever gauge choice would return us to a minimal situation. We 
have investigated two options: 

i) Consider F^M^^'^F^ with a background dependent but non-differential operator M . 
This cannot yield a minimal graviton wave operator because M can only be the Ricci, 
not the Riemann, tensor. 

ii) Consider adding '^R^'"^'^hpa-^v to F^ with new gauge parameter C,. The new terms 
in F"^ then indeed affect the Gauss-Bonnet contributions to the quantum action. However, 
we have verified that no choice of the parameter C will cure the wave operator and make 
it minimal. 

Thus, we are inevitably led to a non-minimal wave operator. This implies that at 
quadratic order in the quantum fields (jf , the action takes the form 

^(2) ^ j dv {^P.j^W^i + \X,,c^''c^ + \Wi:;'D^cl>W,^) (11) 

where P, X and W are symmetric matrices which are in general background field depen- 
dent. Without loss of generality, we may also assume that W^^ ~ W^^ . The covariant 
derivatives satisfy 

[D,,D^]\ct>'' = Y^^\4>'' (12) 

The tensor W is dimensionless and the corresponding vertex can be inserted arbitrarily 
often into a one loop diagram without changing its degree of divergence. This is the 
reason why no simple algorithm which generates all one loop divergences for non-minimal 
operators as in Eq (jlip is known. 

Some partially successful attempts were made though, on which we now comment. 
Pronin and Stepanyantz [24j relied on the diagrammatic methods of 't Hooft to derive 
what they call "master formulas for the divergent part of the one loop effective action 
for arbitrary (both minimal and non-minimal) operators of any order in 4-dimensional 
curved space". However, the authors of [53] assume that W, called K in their work, is 
covariantly constant (see the remarks following Eq (28) of fl^). Basically, they assume 
that the tensor W is always made from products of metric tensors, but this is not the case 
in the present situation. Avramidi [H] and Avramidi and Branson [13 HJ] systematically 
extended the Schwinger-DeWitt method to non-minimal operators, called non-Laplace 
type operators there, and diagonal values of the heat kernel coefficients ao and ai (but 
not yet 02) were giver(|. In particular, the so-called commutative limit defined in Eq (4.2) 

-fSuch a gauge was first considered in |21j in the context of Einstein gravity. There, it was erroneously 
claimed that the choice of C, could affect the numerical coefficient of the well-known non-renormalizable 
two loop on shell divergence of this theory. See 1221 and 1231 for the correct statements. 

§In | 25| . [ai] is proposed as a possible non-abelian generalization of gravity. 
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of [IS] coincides with our approach, but also in 17\ it is assumed that W, caUed a there, 
is covariantly constant, which as aheady noted is not the case in the present situation. 
Already in 1985 Barvinsky and Vilkovisky [13j perturbatively extended the Schwinger- 
DeWitt method to non-minimal operators at one loop order. Their main application was 
to operators of the kind SpD — AV"V/3 for vector fields, where the second term is due to 
non-minimal gauge fixing. Non-minimal operators for gravitation were also covered. We 
have a similar situation, but the offending contribution generated by the Gauss-Bonnet 
term is already of first order in the background curvature. In the parlance of [13j . our 
non-minimal term has background dimensionality 0{l/£'^). This allows us to treat the 
W-terui as a pertubation which simplifies our analysis a lofl For a flat background 
spacetime and to second order in W, an algorithm was first provided Ichinose et al. |14j . 
This work seems to have gone unnoticed so far. No special properties of W were assumed 
in [T3]. We have checked and confirmed this algorithm and extended it to a curved but 
Ricci fiat space plus some terms cubic in W. 

Since we will treat the T4^-term perturbatively, we choose to define the propagator 
G^^{x,x') as the inverse of the minimal part of the wave operator 

{PD' + XUG''^ =-V<5 (13) 

Comparing S^^) + Sfi^, Eq ^ plus with Eq HI]), we read off that for Ricci flat 

background metric we have 

X-V - -^oS^^ + + ^^i^U''\p.) - 2C/-(/„) - 2U'^^^p (14) 

Here, the flrst two terms in W still need to be symmetrized under pair interchange 
fiv <-> pa. In fact, W is then totally symmetric under interchanges kX fiv pa. In 
the second expression we defined 

Note that after contracting a pair of Weyl tensors twice with each other one can always 
put the indices in this standard order [26 . Tensor U has the following symmetries 

Up,iyp(j Upu^u Uiip^up (1^) 

Furthermore, the field equation allows us to require U to be traceless on any pair of its 
indices, i.e. we will drop triple or fully contracted pairs of Weyl tensors. We note that 

{W^-Y\^^{Wp^r\'^ = A{A-d)X^C\;^^ , V^{W\Y\^^Q (17) 

which will be essential in simplifying our analysis. In d = 4, the Gauss-Bonnet term is 
topological, hence it is invariant under any change of the metric and in particular under 
conformal transformations. This explains the vanishing of the traces. The vanishing of 



^Cf the analysis in sect. 3 of 1131 where it is shown that it is better to analyze non-minimal gauge 
theories via the method of Ward identities. This leads namely to powers of the extremal, i.e. the 
background field equation. 
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the divergence of W has its origin in the general coordinate invariance of the Gauss- 
Bonnet term which was preserved by our choice of gauge. 

In principle, we should use the configuration-space metric and its inverse to raise 
and lower pairs of indices on the various tensors. However, in Ricci flat spaces there is 
effectively no difference between the various forms of X (this assumes we drop three-fold 
contractions of two Riemann tensors in the part of X arising from the Gauss-Bonnet 
term, but such a contraction can anyhow be rewritten via the field equation ([5]) and does 
not contribute to the scalars). The same is true for V^^. In d = 4 it is also true 
for W^jy because it is traceless there. In particular, the terms in X involving a pair of 
three-fold contracted Weyl tensors can be dropped. 



4 One loop divergences for Einstein- Gauss-Bonnet 

To find the one loop divergences for the Einstein-Gauss-Bonnet action, we view it as an 
action of the form presented in Eq (|Iip and treat the VF- vertex perturbatively. To zeroth 
order in W and in a Ricci flat space, the divergent part of the one loop effective action 
for a minimal wave operator as in Eq (|I3p is known to be given by 

iTr(lnG),^, = _L_ Tr [a,] 

= Y^Tr(iX^ + ^y2 + ^C2) (18) 

Here e — A—d, [02] signifies the diagonal value of the second heat kernel coefficient [27], Tr 
denotes the functional trace operation, = Y^'^Yf^^, and = C^'^'^'' C^^pa- Actually, 
Eq implies that only the X'^ term contributes to divergences. 

To find the first order effect of the ly-perturbation, we insert one VK-vertex in the 
one loop graph and determine its divergent part in a Ricci flat space, namely 

\l:l{W^^''D^,GD,)^^ = Tr (W'' {D^[D,ai] - [D^D^ai] + ig^.[a2]) 

= Tr {}W^-{D,D^X + '^Y^, ^C^^^^C^P'^r 

+ \W^^ [X^ + \D^X + + .1.^2)) (19) 

Due to Eq (I17p and the remarks at the end of section 3 only the W'^'^Y^ ''^pi, term is 
relevant. 

Inserting two M^- vertices and again finding the divergent part yields 

iTr {W^^'^{x)D,G{x,x')D p>WP'-' {x')D„,G{x' ,x)D ^.^ 

= Tt {^[W^'^'alWP'''' D„,aoD^]Gi^,p,Go + ^[W^'^aoWP'''' D^.a^D p]Go,.p'Gi 
+ [W^''aiWP'^'aoD^]Gi^,p,Go,a' + [W^'aoWP'^' a^D ^]Go^,p,Gi,^. 
+\[W''aoWP'^' ai]Go,,p'Gi^^„, + \[WP'' a^WP a2]Go,,p.G2,^,a' 
+\[W^^''alWP'-'a^]G^^,p>G^^p„,)^.^ (20) 

Here we inserted a heat kernel expansion for each Green function, see [3S], and 
distributed the derivatives, keeping only such terms which in the end can contribute to 
. In particular, we dropped W^'^ , D^WP^ , D^X and D^D^X . A list for the divergent 
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products appearing here can be found in [25] , We repeat here the few ones that are 
actually needed (see Eq (D.19) of [55]) 

Go^^iy'Gi^p — \{gfj.i,V p — 2gi,pV ^ — 2gp^pW^)6 

Gi,^i.'Go,p = ^{2gp^Vp - gi,pVfj, - g^,pVy)5 (21) 

Go^fj,u'G2,pa' = 2Gi,p,^'Gi,pcr' — \g{pijgp(T)5 

Gi.f,i,'Gn^pa' = ^G^(p^)^^ + VV(5 -terms 

where we omitted a factor (167r^e)^^ on the far right hand sides. After substituting 
these expressions, we partially integrate the covariant derivatives off the (5-functions and 
perform the integration over x' . The final result for the relevant terms is 

iTr {^•'D,GDp,WP'^'D,,GD^)^.^ 

1 Tr (^M^^V W^u + T.W''. XW\, X + ^W^, D^W^^ X (22) 



167r^e 

^W^^W^pXYP^ +±WP,XW^pYP^ +^WP,WP,XC^p-'^^ 



where the right hand side is local. The last term explicitly involves the curvature and is 
absent from the flat-space algorithm presented in |T3] . We have verified and agree with 
the complete algorithm of [T3]; cf their Eq (15). Note that the sign of Y in the present 
study is opposite to that in 

Finally, because X contains a term of zeroth order in the Weyl tensor, see Eq (fH]) . 
we also need the two W^X^ invariants, whose coefficients are easily determined. 

Adding the various contributions and keeping only terms which can produce G^, our 
result is 

^^iv = Tr[iX2 + liyM, Y-p YP^ + ^W>^, W\ X^ + l^W>^, XW\ X 

^\W^,W-pXYP^ +±W^^^XW^^YP^ +±WP,D^W^^X 

+ ±WP, WP„ XC^p-- + ^WP, W% WP^ X^ + ^WP, W% XWPp X] (23) 

Calculations with the Mathematica [29] package MathTensor [30] yield in the graviton 
sector 

tr[X2] = -12AiA2(G|,-2G^) 

triW^, Y'^pYP^] = - 12 A2 (G^ - 2G3,) 

tT[W^^W''^X^]^tT[W^,XW''^X] = 576AoAiA2(G^-2G^) 

tr[WP,W%XYP^]=tr[WP,XWpYP^] = -72 X.XUC^ ~ 2Cl,) (24) 

trlWP^D^WpX] ^ 

tr [W'^,WP,X]Gpp'"' = lU\o\l{C^-2Cl,) 

tr [WP^ WP^ X^] = tr [WP, XWP^ X] = 2160 XlXHC^, - 2Cl,) 

Note that each invariant by itself is proportional to Eq. The abbreviations Cy^ and G|,j 
were defined in Eq Q In the here chosen DeWitt-Feynman gauge, the ghost fields do not 
contribute to the on shell G^ divergences. Substitution of Eq ([T4|) into Eq (|23| produces 
our final answer 



(1) ^ ^2 

167r2e 



2 + 6Ai + I2A0A2 - 72A0A1A2 - I8OA0A2 



(G^ - 2G^) (25) 
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This is our main result. It shows that the addition of the Gauss-Bonnet term to the sum 
of the Einstein-Hilbert and cosmological terms induces the on sheh six-dimensional Euler 
density as a new divergence. 

One might conjecture that the divergences would take the form of the dimension- 
ally continued eight-dimensional Euler density Eg- In a Ricci flat space and using the 
equation of motion Eq ([5]), Eg is proportional to the expression graphically represented 
in figure [31 

.(S)-® -.(g) -CS) -cM) 

i^i fii f^i /-~f4 

'^W '-'PI '-'Kl '^K2 

Figure 3: Graphical representation of the Euler density Eg, 

If we limit our attention to the Af sector we find that only the third and fourth terms 
in Eq (|23p contribute. These two terms yield a divergence proportional to 

6C^ + 6Ct,^17CU~7Cj,, + 10Cj,^ . (26) 

This shows that, in contrast to the Eq results, the divergences in the Af A2 sector do 
not take the form of the Euler density Eg. 



5 Conclusions and discussion 

We have demonstrated that the addition of a Gauss-Bonnet term to an action consisting 
of the Einstein-Hilbert term plus cosmological term induces new on shell divergences 
at one loop order of the form of the six-dimensional Euler density Eq. From a strict 
four-dimensional point of view this implies that the Gauss-Bonnet term has no influence 
on the one loop renormalizability of gravity. We thus lay to rest this issue, raised long 
ago by Capper and Kimber [1], at least for the divergent parts of one loop diagrams. 
The question remains open for finite one loop scattering processes and also for a possible 
influence of the Gauss-Bonnet term on the two loop divergences of gravity. Still, given 
the topological nature of the Gauss-Bonnet term in d = 4, our result is as expected. Ex 
nihilo nihil Gt. 

In retrospect we might say that 't Hooft and Veltman [12] showed that if one starts 
from the Einstein-Hilbert term one induces the Gauss-Bonnet term at one loop. Later, 
Christensen and Duff [31 added a cosmological constant to the classical action with again 
the Gauss-Bonnet term being induced. Of course, power counting suffices to predict 
this, because there is only a single scalar quadratic in the Weyl tensor. This is not 
anymore so if one includes the Gauss-Bonnet term in the classical action as we did. 
There then exist two different cubic scalars, but our calculations show that they appear 
only in the combination corresponding to the six-dimensional Euler density. It leads us 
to the conclusion that from the d-dimensional point of view advocated by Capper and 
Kimber [T] , one loop renormalizability of Einstein gravity requires one to extend it to a 
Lovelock theory. However, preliminary calculations show that at quartic order in the Weyl 
tensor it is not the eight-dimensional Euler density Eg which appears, but rather several 
non-topological invariants which vanish in d = 4. It is conceivable that one can arrange 
for such extra terms to produce just Eg by a special choice of coupling constants. That 
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would constitute an interesting constraint on the parameters of the class of all Lovelock 
theories. So far, such constraints were based on stability [HI [2U] or causality [33]; see 
also [31]. We hope to return to this issue in future. 

Our work can also be seen as an interesting setting in which non-minimal wave op- 
erators occur. We have verified the one loop algorithm of [Mj for such operators and 
extended it to curved Ricci flat spaces, both via 't Hooft's non-covariant method and 
by using the covariant Schwinger-DeWitt method. Our result does not depend on the 
precise coefficients in the algorithm: Every contribution individually reduces to the six- 
dimensional Euler density. 

Possibly our work can be continued to all orders in the non- minimal term by extending 
the methods of [HI [TBJ [171 [21] to the case of divergenceless rather than covariantly 
constant W^-tensor. 
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